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Abstract 



Consider an SDE on a foliated manifold whose trajectories lay on compact 
leaves. We investigate the effective behaviour of a small transversal perturba- 
tion of order e. An average principle is shown to hold such that the component 
transversal to the leaves converges to the solution of a deterministic ODE, ac- 
cording to the average of the perturbing vector field with respect to invariant 
measures on the leaves, as e goes to zero. An estimate of the rate of convergence 
. is given. These results generalize the geometrical scope of previous approaches, 

including completely integrable stochastic Hamiltonian system. 
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1 Introduction and set up 

m 

Generally speaking, the original heuristic idea of an averaging principles refers to an 
intertwining of two dynamics where one of them is, in some sense, much slower and 
I/-) ' is affected somehow by the other faster dynamics. An averaging principle in this 

case refers to the possibility of approximate, in some topology, the slow dynamics 
considering only the average action or perturbation which the fast motion induces 
on it. These ideas have appeared long ago and, as mentioned by V. Arnold [3, 
p. 287], they were implicitly contained in the works of Laplace, Lagrange and Gauss 
on celestial mechanics; literature on the matter can be found e.g. among many others 
in [3], Sanders, Verhulst and Murdoch [13J and references therein. Presently, on what 
regards stochastic systems, averaging has been quite an active research field on which 
there is also a vast literature on the topic. Interesting quick historical overviews can 
be found in X.-M. Li [9, p. 806], Kabanov and Pergamenshchikov [7J Appendix], [T3l 
Appendix A]. Among many other works somehow related to the topic, we refer to 
Khasminski and Krylov [5] , Sowers [2] , Namachchvaya and Sowers [TO] Borodin and 
Freidlin [4], [7] and references therein. 

The specific problem that we address in this article is a perturbation of a diffusion 
in a foliated manifold M such that the unperturbed random trajectories lay on the 
leaves. The perturbation are taken transversal to the leaves of the foliation. Here, 
the slow system is the transversal component and the fast system is given by the 
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rescaled y\ , where y\ is the solution of the original SDE perturbed by a vector field 
eK. 

Our results generalize the recent approach by X.M. Li [9] on an averaging princi- 
ple for a completely integrable stochastic Hamiltonian system. In that article, as in 
the classical approach, see e.g. [3], Li has explored the benefits of a well structured 
geometrical coordinates in the state space given by the coordinates of the Liouville 
torus; these benefits include vanishing Ito-Stratonovich correction terms besides also 
vanishing covariant derivative of Hamiltonian vector fields in tangent directions to 
the leaves. We prove here that an averaging principle also holds in a generalized 
geometrical scope, so that this averaging phenomenon occurs independently of sym- 
plectic structures (but with possibly slower rates of convergence). Comparing to Li's 
previous result [9] Lemma 3.2], where the estimates contain a term of order 
our corresponding estimates in Lemma l3.1l are continuous at e = t = 0. Some of the 
rates of convergence of [H Lemma 3.1] is recovered as particular cases in Corollaries 

I2.2l and l2.3l In the main result, we show that in the average, the approximation goes 

_§_ 

to zero faster than | In e| p , for (3 G (0, 1/2). 

The set up. Let M be a smooth Riemannian manifold with an n-dimensional 
smooth foliation, i.e. M is endowed with an integrable regular distribution of di- 
mension n (for definition and further properties of foliated spaces see e.g. the initial 
chapters of Tondeur [15] : Walcak [IB] among others). We denote by L x the leaf of 
the foliation passing through a point x £ M. For simplicity, we shall assume that the 
leaves are compact and that each leaf L x has a tubular neighbourhood U C M where 
U is diffeomorphic to L x x V, where V C R d is an open bounded neighbourhood of 
the origin and d is the codimension of the foliation. We shall assume an SDE in M 
whose solution flow preserves the foliation, i.e. we consider a Stratonovich equation 

r 

dx t = X (x t )dt + Mxt) o dB\ (1) 
i=i 

where the smooth vector fields are foliated in the sense that Xi(x) G T X L X , for 
i = 0,1, ... ,r. Here Bt = (B^, . . . , Bl) is a standard Brownian motion in R r with 
respect to a filtered probability space (O, T%, T, P). For an initial condition xq, the 
trajectories of the solution x% in this case lay on the leaf L XQ a.s.. Moreover, there 
exists a (local) stochastic flow of diffeomorphisms Ft '■ M — > M which restricted to 
the initial leaf is a flow in the compact L xo . 

For a smooth vector field K in M, we shall denote the perturbed system by yf 
which satisfies the SDE 

r 

dyf = X (y e t )dt + Y, My!) ° dB\ + eK{y e t ) dt, (2) 
i=i 

with the same initial condition j/q = xq of the unperturbed system xj. 

Our main result, Theorem 14.11 says that locally the transversal behaviour of 
y% can be approximated in the average by an ordinary differential equation in the 
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transversal space whose coefficients are given by the average of the transversal com- 
ponent of the perturbation K with respect to the invariant measure on the leaves 
for the original dynamics of Equation (pQ). The reader will notice by the end of 
the proofs that compactness of the leaves in fact can be substituted by some other 
boundedness conditions, added also to some rather technical adjustments which we 
will not address here. In the Sections 2 and 3 we present the main lemmas. The 
main result appears in Section 4, where we also present a simple illustrative example. 
In particular, under some symmetry hypothesis on a foliated system embedded in 
an Euclidean space, we use the main theorem to conclude that Lyapunov exponents 
in the transversal direction must tend to zero as e goes to zero, cf. Proposition 14.21 

2 Preliminaries results 

Our coordinate system. Given an initial condition xq G M, let U C M be a bounded 
neighboorhood of xq which is diffeomorphic to L XQ x V and whose closure U C M. By 
compactness of L Xo , there exists a finite number of local foliated coordinate systems 
ipi : Ui — )• Wi x V C R™ x R d , where Wi and V are open sets, say with 1 < i < k 
and xo £ Ui such that: 

1) U = U* =1 U i ; 

2) The leaf L xo = U* =1 ¥> -1 (Wj x {0}), i.e. each Ui is diffeomorphic to the product 

of an open set (with the induced topology) in the leaf L Xo and the vertical 
component V; 

3) If a pair of points p £ U and q G Uj in U belong to the same leaf then their 

transversal coordinates in V are the same; i.e. n((fi(p)) = 7r(^-(p)) where tt is 
the projection on the transversal space V; 

4) For i = l,...,k, ifi has bounded derivatives (obtained reducing open set U if 

necessary) . 

Note that for a fixed y 6 V, the finite union U k ip^ 1 (Wi, y) is the leaf L 
for any x G W%. Natural examples of this scheme of coordinates systems appear if 
we consider compact foliation given by the inverse image of submersions: values in 
the image space provide local coordinates for the vertical space V. 

Next lemma gives information on the order of which the perturbed trajectories 
Vt approaches the unperturbed xt when one varies e and t in equation ([2]); it will 
be used to prove that the dynamics of the rescaled system y\ is such that its time 

average for any function g in M approximates the time average of the spacial average 
of g on the leaves, Lemma 13.11 An exponential factor in the estimates is expected, 
as trivial linear examples show. 

We shall denote the coordinates of a point p G U\ by <pi{p) = (u,v) G R" x H d . 
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Lemma 2.1 Let r e be the first time the process y\ exits the foliated coordinate neigh- 
bourhood Ui as above. For any locally Lipschitz continuous function f : M — > R and 
2 < p < oo we have 



E( sup \m)-f(x.)\ p 

S<t/\T E 



1 
p 

- r v 



where K\ , K<i > are constants depending on upper bounds of the norms of the 
perturbing vector field K, on the Lipschitz coefficients of f and on the derivatives of 
Xq,X\ ■ ■ ■ ,X r with respect to the coordinate system. 

Proof: Initially write xt and y\ , the solutions of Equations (pQ) and (|2|) respectively, 
according to the foliated coordinates ip\. So we write (u t , v t ) := (pi(x t ) and := 

¥>i(yf)- Then 

\f(yt)-f(x t )\ = \fo^(u^vl)-fo^\u u v t )\ (3) 
< C \u\ — u t \ + C \vl — v t \ , 



for some constant C > 0, using the fact that U\ is relatively compact. 

We shall denote ut = (u\, • • • , u"), vt = (v\ , • • • , vf ), u\ = (uf , • • 
(f j' 1 , • • • , v\ ,d ). We also split the horizontal and vertical component of the per- 
turbing vector field K = (K U ,K V ), into coordinates: K u = (if^, . . . , if", ) and 
K v = (Kl...,Ki). 

In our coordinate system, the equations of the horizontal and vertical components 
of the perturbed system u\ and v\ are given by 

T 

duf = Y^bU^^^Bt+boi^v^dt + eKiulv^ (4) 
fc=i 

dv? = eKi(ul,vl)dt, (5) 

with i = 1,2, ■■■,n and j = l,2,---,d, for the induced vector fields bo,bi, . . .b r 
which, together with their derivatives, are also bounded. From equation ([5]) we have 

sup |v| — v s \ < e sup / \K v (u e s , v e s )\ ds 

s<tAr e s<thr e Jo 

< e t sup \K v (x)\ = C x e t, (6) 
xeu 

where C\ = sup^g^ |K(x)|. From equation (J4j) we have in each i-th component, for 
s < T e : 



u e f - ui 



= E f ( b U U >r) - &U«r,Vr)) ° ^ r fc + (7) 

/ - b l (u r ,v r ))dr + e / K l u (u e r ,v e r )dr. (8) 

Jo io 
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In terms of ltd integral, 



Jo 



(&l«><) -bi(u r ,v r )) odB. 



(&fc(«r>«r) -bi(u r ,V r ))dB, 



+ o / I V6 1 • b k(u € r ,v e r ) - W fe ■ b k (u r ,v r )]dr. 



Hence, taking the absolute values in both sides of Equation (jSJ) we get, for each i: 



\u e ' v -ul\ < 



E 

k=l 

-t 

fc=l 



ib\.{u e r ,vl) — b\(u r ,v r ))dB r 



+ 



\Vb\ ■ b k {u e n v e r ) - Vb\ ■ b k (u r ,v r )\ dr 



o 



+ / \b (u e r ,v^.) - b (u r ,v r )\dr + e \K l u (ul,v € r )\dr. 
'o Jo 



(9) 



Functions b Q and (V6^ • b k ) are Lipschitz, hence for a common constant C2, 



f\mif r y r )-v k {u r ,v r ))dB k r + 
k=l J ° 

C2 / \v* — v r \ dr + C2 / \u e r — u r \ dr + es svly>\K u \ . 
Jo Jo u 



(10) 



The first deterministic integral, together with inequality (0) yields: 



<•'-< < 



E 



+C2 / |u£ — u r \ dr + Cies. 
Jo 

Now, for p > 1, there exists a constant C3 such that 



Sjj] [ S (bi(u<r,V* r )-bi(u r ,V r ))dB 
k=l J ° 

+ C 3 Ci(^K-u r \dr\ +C 3 (C 1 es) p . 
Cauchy-Schwartz inequality yields: 



2^P 



+ C 3 (C^es 2 ) 



CsJ] l\bl{u%y r )-bl{u r ,v r ))dB 

+C ?J Cl sP' 1 [ |< - u r | p dr + C 3 (des 
Jo 



+ C 3 (CiC7 2 es 



Hence, 



E sup \u e /-ui\ P < C 3 E sup V 

s<Mt< s<tAr e , -, 

— — fc=l 



(b{(u e r ,v e r ) -bi(u r ,v r ))dBl 



+ C 3 (dC 2 et 2 ) 



2\P 



+C 3 C P t p ~ l E sup ( |< - u r | p dr + C 3 (Ciet) f 

S<t/\T e JO 



b l k (ur,v r )) 2 dr 



k=l u ° 

+C 3 C%i p - 1 [ E( sup |<-u r | p ) dr + C 3 (Ciet) p 

JO \s<rAT ( / 



2\P 



C 3 (CiC 2 ei 2 ) 



where we have used classical L p -inequality for martingales (e.g. Revuz and Yor [11]). 
Using again the Lipchitz property of each b k for the terms in the brackets above: 



r f(AT £ 

E / 



uj:) - b l k (u r ,v r )) 2 dr 



< 2Cl 



tf\T e 



tf\T e 



v% — v r \ dr + / |u^, — u r | dr 



o 



< 2C| ( / CfeV + 



iAr e 



sup | u\ — u r I dr 

s<rAr £ 



< CiCfeH 3 + 2Cf / sup 

JO s<rAr £ 



llif; — u r \ dr. 



(11) 



We end up with: 



rt/\T 

E sup \uf-v? s f < C 4 VE C 2 C' 1 W + 2C 2 / sup 

s<Mt ! ' , =1 L JO s<rAi 

+C 3 C%t p - 1 f if sup K-n r | p )dr 

J \s<rAr e J 



u% — u T I dr 



p/2 



+C 3 (C 1 C 2 ei 2 ) p + C 3 (C ie t) p . 



(12) 



For p > 2 one can use Cauchy-Schwartz again to conclude that there exists a constant 
C5 such that the last expression is less than or equal 

C 2 C 1 et 3 / 2 ) +C 5 CftV / E sup \ul-u r \ p dr + C 3 {CiC 2 et 2 ) p 

' Jo s<rf\T e 

+C 3 C P tP- 1 !\( sup |< - U A P ) dr + C 3 (detf 

Jo \s<rAr e J 



C 5 (ddet^Y + C 3 (CiC 2 et 2 ) P + C 3 (detf 



+ (c 5 c p t E ^ + c 3 c p t p - 1 



E I sup \u e r — u r \ p ) dr. 

\s<rAT e 



Now, summing up over i in the inequalities above leads to 



E sup \u e r -u r \ p < C 5 (c 2 Ciet 3 / 2 ) P + Cz(C 1 C 2 et 2 ) p + C 3 {C 1 et) p 



!<(AT< 



+ (c 5 c p t^ + c 3 c\t p ~ x \ l\( sup 

^ ' Jo \s<rATt 



\u e r — u r \ p } dr. 



We use now the integral form of Gronwall's inequality to find that: 

E ( sup |< - UrA < C 6 e p t p (1 + t p ) exp{C 7 (t p/2 + t p )} 

\s<iAr £ / 

< C 8 e p t p (l + t p )exp{C 9 t p }. 
Going back to the inequality now we have 

\f(yi) - f(x t )\ p < Cio K - + Cio K - u t | p 

hence: 

e( sup 1/(^-/(^)1^ < CiqE sup |«J-« a |P + Ci E sup |<-u s | p 

\s<tAr E / s<tAr E s<(At< 

< Ciie p t p + Ci 2 e p t p (l + t p )exp(C7 9 t p ) 



< Ci 3 e p ^(l + t p )exp(C 9 
From here, finally, one concludes that there exist constants K\ and if 2 such that 



E I sup \f(yl)-f( Xs )\n <K ie te W (K 2 t p ). 

Vs<tAr E / 

□ 

Next corollary includes the case of a completely integrable stochastic Hamiltonian 
system when one uses the action-angle coordinates, cf. X.-M. Li p2 Lemma 3.1]). 

Corollary 2.2 If the vector fields Xq, ■ ■ ■ ,X r depend only on the vertical coordinate 
(null derivative in the directions of the leaves, as in the Hamiltonian case JE/) then 
the estimates above can be improved, and for p > 1 there exists a constant K\ such 
that 

e( sup \f(y e s ) - f(x s )\ p )] P ^K^it + t 2 ). 

\s<tAr e / . 
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Proof: In this case the correction term of the Stratonovich stochastic integral in 
terms of Ito integral in inequality ([9]) vanishes, and also so does the determinist 
integration of \u e r — u r \ in inequalities (jlOp and (jlip . Hence inequality (|12p improves 
to 

E sup \uf-u\\ p < cJ^det^Y + C 3 (C 1 C 2 et 2 ) p + C 3 (C ie t) p . (13) 

The argument in the rest of the proof follows straightforward for p > 1 skipping 
Gronwall inequality. 

□ 

Next corollary includes the case Xq = 0, cf. [9j Lemma 3.1(2)] for stochastic 
Hamiltonian systems with action-angle coordinate system. 

Corollary 2.3 If in addition to conditions of Corollarv \2.2\ above, we have that the 
deterministic vector field Xq is constant when represented with respect to a certain 
local coordinate system in Ui (i.e. bo has null derivative w.r.t u and v) then, for 

3 

p > 1 the estimates can be improved further to K\e{t + ). 

Proof: Besides the vanishing terms already mentioned above, the second determin- 
istic integral on the right hand side of inequality ([9]) also vanishes. Hence inequality 
(|12p simplifies further to 

E sup \u e /-ui\ p < C 5 (CiC^^Y + C z {Cietf . 

s<iAr' ^ ' 

□ 

Yet, from the proof of the Lemma 12. II we have the following 

Remark 2.4 For 1 < p < 2 and t sufficiently small, there exist constants K\ and 
K2 such that 

< Kiet exp(K 2 t p ). (14) 

Proof: One can no longer use Cauchy-Schwartz after inequality f)12[> . Alternatively, 
from (I12p . use that 

E sup \u e s - u s \ p < C 5 (c 2 Ciet 3 / 2 ) P + C 5 CfW 2 E sup \uf - u\\ P + C 3 {C l C 2 et 2 ) p 

s<Mr< ^ ' s<1At' 

+C 3 Cf / e( sup lu^-UrA dr + C 3 C p e p t p . 

JO \s<rAr e J 

If we fix an < 5 < 1 and take t sufficiently small such that 1 — C^C\t p l 2 > 5 then 
E sup \ul-u s \ p < 6- 1 C 5 (c 2 C 1 et 3 / 2 Y + S- 1 C 3 {C 1 C 2 et 2 ) p 

s<1At< ^ ' 

+5- 1 C 3 C p t p - 1 f\( sup |< - u r \ p \ dr + 5- l C 3 C p e p t p . 

JO \s<rAr £ / 



E| sup \f(y e s )-f(x s )\ p 

S<tAT e 
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And one completes the calculation as before using the integral version of Gronwall 
inequality. 

□ 

3 Averaging functions on the leaves 

Consider a differentiable function g : M — > R. The leaf L p passing through a point 
p £ M contains the support of an invariant measure \x v for the unperturbed system 
([T]); we shall assume that [i p is ergodic. We shall work with the following function 
defined for each leaf Q 9 : V C H d — > R given by the average of g with respect 
to these measures on the leaves. Namely, if v is the vertical coordinate of p, i.e. 
p = ip(u, v), then: 

Q9( v ) = I g(x)dn p {x). 

J Lp 

We assume that Q 9 has some degree of continuity with respect to v. This assumption 
appears in two levels in Lemma 13.41 1) Riemann integrability of Q 9 {Ti{y^)) with 
respect to r guarantees the convergence to zero; 2) a-H61der continuity guarantees 
the rate of convergence. 

Our next step is to estimate the time average of g and Q 9 along the perturbed 
system y\. To use the notations and results we have introduced before in local 
coordinates, we shall write simply ir(p) = v for the composition of the projection on 
the second coordinates with the local chart ip(p) = (u,v). 

Here the stopping time r e denotes the first exit time of the open neighbourhood 
U C M which is diffeomorphic to L XQ x V. We have the following estimates for the 
difference of the averages of functions g and Q 9 . 

Lemma 3.1 Given a function g : M — >• R let Q 9 : V — s- R be its average on the 
leaves. For s,t>0 write 

/•(■S + t)AET e 

5(e,t)= / gm-Q g Wr_)) dr. 

Then 5(e, t) goes to zero when t or e tend to zero. 

Moreover, if Q 9 is a-Holder continuous with a > then for p > 1 and any 
P £ (0, 1/2) we have the following estimates: 

^Esup|<5(e,s)| p y < Vt\\ne\~vh(t,e), 

where h(t, e) is continuous for t, e > and converges to zero when (t, e) — > 0. 
Proof: 

The proof consists of considering a convenient partition of the interval (s/e A 
t € , (s + t)/e A t € ) where we can get the estimates by comparing in each subinterval 
the average of the flow of the original system (on the corresponding leaf) with the 
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average of the perturbed flow (possibly transversal to the leaves). These estimates 
in each subinterval are obtained using Lemma 12.11 So, a key point in the proof is a 
careful choice of the increments of such a convenient partition. For sufficiently small 
e, we take the following assignment of increments: 



At 



(s + t) A r £ - s A r e 
hie p 



Hence, the partition t r , 



A r e + nAt, for 1 < n < N — 1, is such that 



s s + t , 

- A T € = t < *1 < • • • < tN-1 < A T e . 

e e 



with N = N(e) = [e~ \ lne| p ] where here [x] denotes the integer part of x. 
Initially we represent the left hand side as the sum: 



n=0 



-Ar e 



N-l r t n+1 

el g(y e r )dr = e^2 M)dr + e I g{y e r )dr. 



Denote by 9t the canonical shift operator on the probability space. Let Ft(-,u>) 
with t > be the flow of the original unperturbed system in M. Triangular inequality 
splits our calculation into four parts 



where 



|<5(e,t)| < |Ai| + |^ 2 | + |A 3 | + |A4| 



(15) 



A, 



A A 



€ 

n=0 
N-l 

( 

n=0 
N-l 



tn+l 



E 

n=0 
N-l 

E 

n=0 
-1 

^2 eAtQ9(7r(yl)) 



[g{yl)-g{F r . tn (yt,6t n ^)))\ dr, 



g(F^ tn (yl,e tn (cj)))dr - AtQ<> (ntfj) 

(s+t)Aer e 



n=0 



Q 9 (n(y%))dr, 



sAeT e 



-Ar e 



g(Vr)dr. 



We proceed showing that each of the processes A\,A2,Aa and A4 above tends to 
zero uniformly on compact intervals. In what follows, we will explore many times the 
fact that for a > and b G R, e a \ In e\ b goes to zero as e \ 0. Hence, by construction, 
except when T e < s (where, restricted to which the lemma is trivial) both At and iV 
go to infinity when e tends to zero. 
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Lemma 3.2 Process A\ converges to zero uniformly on compact intervals when e 
goes to zero. More precisely, we have the following estimates on the rate of conver- 
gence: For any 7 6 (0, 1), there exists a function h\ such that 



Esup|Ai| p < K x te^ h x {t,e) 

s<t ) 

where h\ is continuous int,e> and converges to zero when (t,e) \-t (0,0). 

Proof: Initially note that by triangular inequality, and putting the supremum inside 
the integral we get 



E sup I Ai| 



s<t 



N—l / r nt 

< «E E / su p Ws)-g{Fs-t n {yt,o u M))\dr 

n=0 V Ut„ t n <s<r 



If - + - = 1, by Holder inequality we have that the estimate above is again bounded 
by P 



N-l 



E ■ 



71=0 



tn+1 



dr 



tn+l 



sup \g(yt)-g(F s _ tn (yl,e tn (u)))\ p dr 

t n <s<r 



2ti 



(16) 

The increment At is a random variable bounded by t\ lne|~. In the estimates 
below At will denote this very deterministic upper bounded. Analogously to N 

which has order [e | In e | »> J. Hence, with this notation, last inequality is again 
bounded by 



N-l 



n=0 
N-l / 



At sup \g(.yl)-g(F s - tn (y e tn ,e tn (uj)))\ 

t n <S<t n + \ 



< eAt ^2 E 



n=0 



sup \g(yl) - g(F s _ tn (y e tn ,6 tn (u>)))\ 

t n <s<t n +\ 



Lemma l2. II and its corollaries says that for all < n < N — 1 above, the function 
g evaluated along trajectories of the perturbed system compared with g evaluated 
along the unperturbed trajectories, both starting at y\ satisfies: 



E sup \g(yl)-g(F^ tn (yl,e tn (u;))\ p 

t n <S<t n +l 



KK^Ate^^ 



Note that here, possibly Lemma [2.11 might be applied using different (finitely many) 
local coordinate systems ipi for each n. 



Esup|Ai| p 

s<t 



< K!Ne 2 (At) 2 e K ^ P 
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^[e^llneff] t 2 \ lne\^e K ^ lne l P > P 
Ktte' hi(e,t) 



4 ' 3 



for any 7 G (0, 1) where 



hi(e,t) =t eV |lne|^ exp \ ( ) In e + K 2 t p | In e |2/; 



which satisfies the required properties for /3 E (0, 1/2). 



□ 



Lemma 3.3 Process A<i in equation [15\) goes to zero with the following rate of 
convergence: 



Esup|^ 2 | p 

s<t 



< Ks/i\ lne| 



for a positive constant K. 
Proof: We have 





1 




Esup|^ 2 | p 


V 




< e 


E 


s<t 







N-l 



n=0 



tn+l 



E [ / 9{F r -tM n A»))dr - AtQ9(n(yl))} 



N-l 



E 



n=0 
AT-1 

= 6 At E 



tn+l 



n=0 



E 



s(F r _ t Jyf ,0 t >)))dr - AtQ°(M)) 



1 ftn + l P 

' g(F^ tn (yl,0M))dr-Q 9 (7r(yl)) 



At 



For all rt = 0, . . . , N— 1, by construction the two terms inside the modulus converges 
to each other when At goes to infinity. Moreover, as in (9, Lemma 3.2] by Markovian 
property and central limit theorem, the rate of convergence has order ^= when At 
goes to infinity. Hence, for small e we have 



Esup|^ 2 | p 

s<t 



< Ke N(At) 



< Ke 



e he p 



1 — 
t\ In e| p 

\/t| lne| p 



KVi\ ln ( 



□ 
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Lemma 3.4 A3 converges to zero when t or e go to 0. Moreover, if Q 9 is a-Holder 
continuous with < 6 < a then the rate of convergence is given by 



Esup\A 3 \ p ) < Ket L+ 2 I In e| ? 

s<t 



for a positive constant K. 
Proof: By definition 



N-l 



n=0 



> / Q 9 W r ))dr 

-Ar e 



and the convergence to zero here corresponds to the existence of the Riemann inte- 
gral. Moreover, assuming that Q 9 is a-H61der continuous for 5 < a; since y\ is also 
a-H61der continuous for a < \ (classical result, see e.g. [H]), then we have that the 
previous expression is again bounded by 



JV-l 



< e^At sup \Q 9 ((7r(yl)))-Q 9 ((7r(yl 

< Ke(At)^N 

< Ket 1+ 2 he p. 



□ 



Lemma 3.5 Process A4 converges to zero with 



Proof: Denoting 



Esup|A 4 | p )' < Cte\\ne\p . 

s<t 



C = sup|e/(x)|. 



The result follows straightforward since 



g(y e r )dr 



< CeAt = Cte\ he p. 



Now, going back to the proof of Lemma [3. 11 Note that each of the four estimates of 

Lemmas l3.2H3.5l allows a factorization which has a common factor lne|~p times 
a continuous function which goes to zero when (t, e) — > (indeed, in Lemma [3T3l use 
a f3' 6 ((3,1/2)). Lemma 13. II now follows by inequality (|15() . 

□ 
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4 An averaging principle 

We state the averaging principle in the next theorem. To use Lemma 13.11 of the 
previous section we have to assume regularity in the average function Q 9 , which 
naturally depends on g, on the foliated coordinate system and on the transversal 
behaviour of the invariant measures on the leaves of the original foliated system. We 
are going to assume the following: 

Hypothesis (H): For any Lipschitz continuous function g on M, its corresponding 
average function Q 9 on the transversal space V which indexes the leaves is also 
Lipschitz. 

This hypothesis holds naturally if the invariant measures \x v for the unperturbed 
foliated system has sort of weakly continuity on p. For deterministic systems it 
corresponds to a certain regularity in the sense that there is no bifurcation with 
respect to the vertical parameter v £ V. 

We use the derivative of each component of 

7r(-) = (7ri (•),..., Tr d (-))ev cn d . 

to get the averages Q dni ^ K \x) of the real functions dni(K), i = 1, . . . , d on each leaf 
L x - 

Theorem 4.1 Assume that the unperturbed foliated system (0]j on M satisfies hy- 
pothesis (H) above. Letv(t) be the solution of the deterministic ODE in the transver- 
sal component V C R" 

§ = (Q* lW ,...,0* r ' w )Kt)) (17) 

with initial condition v(Q) = tt(xo) = 0. Let To be the time that v(t) reaches the 
boundary of V , then 

(1) For all < t < T , (3 £ (0, 1/2) and 2 < p < oo (I < p for small t, cf. Remark 



E (sup 






n (y e S AT* ) - v(s) 


')] 


\s<t 







1 

V r- P 



<Vi\lne\ rh(t,e), 

where h(t,e) is continuous and converges to zero when e or t goes to 0. 
(2) For 7 < 0, let 

T y = M {dtat(i;(t),0y)<7} 
The exit times of the two systems satisfies the estimates 
P(er e < T 7 ) < 7 - p t? | In e\-^h(t, ef . 
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Proof: The gradient of each real function 7Tj is orthogonal to the leaves, hence by 
Ito formula, for i = 1,2, . . . ,d we have that 



d^(K)(y%) ds. 



Lemma |3, II for the function dni(K) in M, triangular inequality and hypothesis (H) 
imply that 



Vi(t) 



< 



tAT E 



Qd MK) ^(yt^.)) -Q d ^ K \v(s))\ ds + \5i{e,t)\ 



< a 









f 




-v(s) 


>0 







where each is the Lipschitz constant of Q d7r d K ) anc [ Si(e, t) is defined in Lemma 
13. 1L Summing up the i's and using Gronwall lemma we have, for a constant C: 

n 

Tr(yl^e) -v(t)\ <e Ct J2\8i(e,t)\. 

i=l 

The first part of the theorem follows by Lemma 13. 11 For the second part we have 
the following estimates 



(er e < T 7 ) < 



sup 

s<T~,f\eT e 



v(s) - 7T (y\ AT * ) 



>7 



< 7 _P E sup 

V s<T-,Aer e 



v(s) - ir (y e tAT ej 



< 



7 



t| lne| vh{t,e) 



< j- p 1 2 \lne\- 15 h(t,e) p . 



□ 



4.1 A detailed example: 

The following simple example illustrates the framework where the averaging principle 
described in this section holds. Consider M = R 3 — {(0,0, z), z G R} with the 1- 
dimension horizontal circular foliation of M where the leaf passing through a point 
p = (x,y,z) is given by the circle L p = {(\/ x 2 + y 2 cos 6, \J x 2 + y 2 sm.6, z), 9 6 
[0, 2ir]. Consider the foliated linear SDE on M consisting of random rotations: 

(0 -1 0\ 

dx t = 1 \x t (Ai dt + A 2 dB t ). 
V / 

For an initial condition po = (xq, yo, Zo), say with xq > consider the local 
foliated coordinates in the neighbourhood U = R 3 \ {(x, 0, z); x < 0; z G R} given 
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by cylindrical coordinates. Hence, using the same notation as before ip = (u, v) will 
be denned by <p : U C M — >■ (— n, tt) x R>o x R, where u G (— vr,7r) is angular 
and v = (r,z) G R>o x R such that : (u,v) h-> (rcosn, rsinu, 2) G M. In this 
coordinates system, the transversal projections tti and ir 2 correspond to the radial 
r-component and the z-coordinate, respectively. 

For Ai,A2 £ R with |Ai| + | A2 1 > 0, the invariant measures [i v in the leaves L p 
passing through points p G M are given by normalized Lebesgue measures in L p , 
which here corresponds to the normalized angle 1-form. Note that Hypothesis (H) is 
satisfied. We investigate the effective behaviour of a small transversal perturbation 
of order e: 

/0 -1 0\ 

dy\ = Ai 1 y\ [\ x dt + A 2 dB t ) + eK(y e t ) dt. 
\ / 

with initial condition xq = (1,0,0). In this example we shall consider two classes of 
perturbing vector field K. 



Constant perturbation. (A) Assume that the perturbation is given by a vector 
field which is constant K = (^1,^2,^3) with respect to Euclidean coordinates in 
M. Initially, to fix the ideas, assume that = 0. Then, not only the average on 
the z-component vanishes, i.e. Q d7T2K = 0, but also, by the geometrical symmetry 
of K with respect to the invariant measure, the average radial r-component also 
vanishes, i.e. Q dniK = 0. Hence the transversal component in the main Theorem 
IO is constant v(t) = (r(0),z(0)) for all t > 0. 

Theorem 14.11 establishes a rate of convergence to zero of the difference between 
the initial radius r(0) = 1 and r(^— -) = TTi(y e tAT (), the radial component of the 

perturbed systems, precisely, that 



E 



sup 

s<t 



S A T e 



1 



^oes to zero as e and t goes to zero with the prescribed rate of convergence. 
For comparison, indeed, the perturbed systems actualy has solution 



/ 



n 



cos ( ^ + X 2 Bt 



\ 



sin ( ^ + X 2 Bt 



\ 







J 



h sin ( M + \ 2 Bt_ ) + k 2 cos ( ^ + X 2 Bt 



Ait 



V 



-k x cos ( ^ + \ 2 Bt_ ) + k 2 sin ( ^ + \ 2 Bt_ ) - ki 







One can easily get an explicit description of the r-component: By normalization 
and using the symmetry, one can fix any k\ and k 2 ; for simplicity, we shall fix 
K = (1, 0, 0), hence, in this case, for t < r e 



r(t) = 1+e 2 



2 + cos ( — + \ 2 Bt_ 

e e 



cos 2 f ^ + \ 2 Bt_ j +sin ( ^ + X 2 Bt_ 
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Now, using that r(t) + 1 is bounded and that \r(t) — 1| = \r 2 (t) — l\/\r(t) + 1| one 
finally finds that 

\r(t) - 1| < K(t)V~e 

for any t < r e , according to the boundaries on the the rate of convergence stated in 
the theorem. 

(B) Vertical perturbations. Now, for constant and vertical K = (0,0, £3), the 
radial average Q dniK is null but Q dn2K equals ks for every leaf in M. Hence the 
averaged system v(t) = (r(0), k^t) is constant in the radial component and increases 
linearly in the z-coordinate. The perturbed systems has the simple solution 

\ 



n 



( cos^ + A 2 Bi 
sin(^ + A 2 £i) 



V 



kit 



I 



Hence, the comparison 



7T2 \ VtKT*- 



v(t)\ = 



for all t > and the convergence of the theorem is trivially verified. 



Linear perturbation. Consider a linear perturbation of the form K(x, y, z) = 
(x, 0, 0). In this case, again, the z-coordinate average vanishes trivially. For the 
radial component, we have that dft\K = ro cos 2 u, where u is the angular coordinate 
of p whose distance to the z-axis (r-coordinate) is ro. Hence the average with respect 
to the invariant measure on the leaves is given by Q dniK = r/2 for leaves with radius 
r. The transversal system stated in the Theorem is then v(t) = (e5r(0), z(0)). Hence 
the result guarantees that the radial part of y \ AT € must have a behaviour close to 

the exponential ea in the sense that 



E (sup 










\s<t 







goes to zero when e goes to zero. 

The fundamental solution of the linear perturbed Stratonovich systems yf is given 
by the exponential of the matrix for each fixed t 

et -(X 1 t + X 2 B t ) \ 

X 1 t + X 2 B t 

/ 

where the eigenvalues corresponding to the first two coordinates (horizontal plane) 
are 

-et±^eH 2 -A{X l t + X 2 B t y 
Ai 9 — • 
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whose real part is given by et/2 with probability increasing to 1 as e goes to zero. 
This exponential rates coincides with that one above guaranteed by the Theorem 

Lyapunov exponents. For foliated manifolds embedded in H N , the symmetry 
of the perturbing vector field K with respect to the geometry of the leaves, hence 
also with respect to Lebesgue invariant measure, as presented in the case of constant 
K, has implied that the transversal average Q dnK vanishes. This phenomenon also 
appears in a couple of other examples where the leaves are not only diffeomorphic 
to each other, but also has this symmetry in the sense that the integration of a 
constant perturbation K 6 R^ with respect to the Lebesgue measure is zero. To 
mention a couple of simple examples: the spherical foliation of R n \ {0}, nested 
torus (increasing the smaller radius) foliation of the solid torus minus the central 
circle S 1 x D 2 \ S 1 C R 3 or more generally (when they exists) tubular foliation of 
R n \ {C}, with C a compact set (this context also includes the Hamiltonian case 
with the Lyouville foliation of the symplectic structure of R 2n as in [9]). In these 
symmetric geometrical configuration, if the invariant measure on the leaves are the 
Lebesgue measures (taking gradient Brownian motion on the leaves for instance, 
as in [6]) the averages of Q dn vanishes, hence our main theorem says that on the 
average, the trajectories of the perturbed system stay somehow close to the initial 
leaf. 

Lyapunov exponent of the system in the direction of a tangent vector v £ T XQ M 
contains the long time behaviour of points close to xo in the direction of v, for details 
on the definition, properties, existence conditions, multiplicative ergodic theory, etc 
see e.g. among many others L. Arnold [2] and the references therein. In particular, 
under the symmetric geometrical circumstances above, if there exists the Lyapunov 
exponents of the perturbed system y\, Theorem 14.11 will imply that in transversal 
directions the Lyapunov exponent can not be too far from zero. This vanishing 
property must happen with multiplicity given by the codimension of the foliation, as 
in the examples of the paragraphs above, where the asymptotic relevant parameters 
(rotation number ([H], p], [5] and Lyapunov exponents) do exist. In short: 

Proposition 4.2 (Continuity of Lyapunov exponents) Assume that the per- 
turbed system y\ does have Lyapunov exponents a.s. at the assigned initial condition. 
If the averaged perturbation in the leaves vanishes, i.e. Q = 0, then a number, 
given by the codimension of the foliation, of Lyapunov exponents in the spectrum 
goes to zero as e goes to zero. 

Proof: In fact, Theorem l4.1l savs that the perturbed system increases in the transver- 
sal coordinates with order |vr(y|)| ~ \fe~t h(et,e). Hence, any exponential behaviour, 
if it exists, must be in the direction of the leaf L XQ . 

□ 
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